Scaling properties of one-dimensional Anderson models in an electric field: 

Exponential vs. factorial localization 
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We investigate the scaling properties of eigenstates of a one- dimensional (ID) Anderson model in 
the presence of a constant electric field. The states show a transition from exponential to factorial 
localization. For infinite systems this transition can be described by a simple scaling law based on a 
single parameter X x = loo/ld, the ratio between the Anderson localization length and the Stark 
localization length i e i- For finite samples, however, the system size N enters the problem as a third 
parameter. In that case the global structure of eigenstates is uniquely determined by two scaling 
parameters A at = lao/N and X^ = ioo/'ei- 
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I. INTRODUCTION 



In recent years several studies have investigated the 
influence of constant uniform electric fields on the local- 
ization of electrons in one-dimensional (ID) systems with 
on-site randomness. In the absence of dc-fields it is by 
now well known that even small amounts of disorder lead 
to an exponential localization of all eigenstatesatl. On the 
other hand, in the case of a single-orbital tight-binding 
model of an electron in a periodic potential, application 
of a static electric field is known to generate a discrete, 
uniformly spaced eigenvalue spectrumtl, known asJStark 
ladder, with all eigenfunction localized factoriallya. For 
weak fields the wavefunctions may be extended over sev- 
eral lattice periods, but with increasing field strength the 
electron tends to be localized on a specific site. This is 
known as Stark localization_and has been observed ex- 
perimentally in superlatticesa which are commonly used 
for such measurementsa. 

In infinite samples the localization of eigenstates can 
be characterized in terms of the localization length; the 
latter is commonly defined from the amplitude decay of 
eigenstates in the limit \n\ — > oo, where n labels the site 
in the tight-binding picture. The most powerful and in- 
formative method available for such studies is the transfer 
matrix method. In the presence of a strong electric field, 
however, it appears to be less efficient due to the facto- 
rial nature of the Stark localization. Moreover, for finite 
systems the structure of eigenvectors cannot be charac- 
terized in the same way. One then needs to use other 
quantities (such as the inverse participation ratio), that 
are valid both for finite and infinite samples. Through 
the use of scaling conjectures, one can link then the prop- 
erties of eigenstates in infinite samples to those of finite 
samples. Since the scaling approach proved to be ex- 
tremely useful in-describing conductance and its fluctu- 
ations (see, e.g. ,013) in the theory of disordered solids, 
it seems natural to use this approach also in order to 
describe localization properties of eigenfunctions of ID 



disordered systems in the presence of constant electric 
field. 

In this paper we study the ID Anderson model in the 
presence of a constant electric field in view of scaling 
properties of its eigenstates. The main question that we 
want to answer is whether the up to now known cqi ' 
alence between quasi-lD and ID disordered modelsE 
can be extended in order to include also systems with 
constant electric field. For this purpose we analyze the 
scaling properties of information lengths for infinite and 
finite samples, which were used successfully ituthe studies 
of one and quasi-one dimensional systemsoLJ. Contrary 
to the standard Anderson case, where the ratio of the 
Anderson localization length and the sample size N, 
i.e. Xn = loo/N, is the only relevant scaling parame- 
ter, we find in the present case an additional scaling pa- 
rameter Aoo = loo /I A- Here l c \ is the Stark localization 
length, which arises from the applied constant electric 
field. Hence, the structure of eigenvectors in our model 
is characterized by two scaling parameters Xn , ^oo ■ 

The structure of the paper is as follows. In Section 2 
we describe the mathematical model and briefly summa- 
rize the known facts about the two limiting cases that 
appear for our model. In Section 3 we discuss different 
definitions of localization length, which are used in our 
numerical simulations. In Section 4 we present numerical 
data on scaling of localization lengths of eigenstates in in- 
finite and finite systems. Finally, in Section 5 we study 
the scaling of the whole distribution of eigenvectors. Our 
conclusions are summarized in Section 6. 



II. THE MODEL 



Our starting point is the ID Schrodinger equation in 
the tight-binding approximation 
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.dlp n (t) . 



neF)^ n (t) + Vil> n+1 (t) + Vi/> n -i(t) , 

(1) 



where ip n (t) denotes the probability amplitude for an 
electron to be at site n. Moreover, e n is the local site 
energy, V is the hopping element, e is the electron charge 
and F the strength of the applied dc field. By applying 
the transformation ^> n (t) = exp(— iEt) (p n one obtains 
the stationary equation 



Eip n = Vifi n+ i + (e„ + neF)if> n + Vifi 7 , 



(2) 



for the eigenvalues E and the corresponding eigen- 
states tp n (E). We can distinguish two limiting cases 
which are relevant for our study: (a) perfect system (i.e. 
e n = e) with non-zero electric field F ^ 0, and (b) zero 
field (i.e. F = 0) with random on-site potential. 
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FIG. 1. Two representative eigenfunctions of the ID tight- 
binding model (|l|). (a) Stark regime with factorial localization 
(W = 0, eF = 0.5, Z e i ~ 2). (b) Anderson regime with exponen- 
tial localization (W = 5, eF = 0, loo ~ 3.4). The dashed line has 
slope 2/lao. 

In the case of a perfect system with F / 0, one can 
prove that the corresponding eigenstates <fi(E), known 
as AjjZ^nnier-Stark states, show a generic factorial decay 



<p n (E) = J m _„(2/eF) - (l/eF)H/(|n|!) 



±oo (3) 



where J n is the Bessel function of order n. Wannier- 
StarkpStates constitute a complete set of energy eigen- 
statesEJ. Their eigenvalues, r-E m — meF form the so- 
called Wannier-Stark laddera. A particular Wannier- 
Stark state <p n is factorially localized around the n— th 
site, with a localization length of the order of l c \ ~ 1/eF, 
i.e. the electric field appears in the denominator of the 
localization length in Eq. (||). This underscores the fact 
that F cannot be treated as a small perturbation to the 
field-free Hamiltonian. An example of such a state is 
presented in Fig. la. 



The other limit of interest corresponds to zero elec- 
tric field with e„ random and (^-correlated, chosen from 
a distribution V e with mean zero and variance a 2 . Be- 
low, in our numerical investigation we will always assume 
that V e is a uniform distribution in [-W/2, W/2]. Such a. 
model is known in the literature as the Anderson modelEJ 
and has been studied in great detail in the context of 
disordered materials. It was shown with mathematical 
rigor that in the limit of infinite samples this model dis- 
plays exponentially localized eigenfunctions, no matter 
how small the disorder is (see Fig. lb). The rate of de- 
cay is measured by the Lyapunov exponent 7 which may 
be evaluated by the transfer matrix method. To this end, 
one has to study the asymptotic behavior of the random 
matrix product Y[ M n , where M n is defined through the 
relation 



= M n £„; M n = 



v n -1 
1 



E-e ri 
V 



(4) 



for the vector = {x n , x n -i) with the matrix M n known 
as the transfer matrix. The localization length Zoo is 
hence the inverse Lyapunov exponent 7; the latter is 
evaluated as the exponential rate of increase of an ini- 
tial vector £1, 
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(5) 



Although the Lyapunov exponent 7 for finite N depends 
on the particular realization, of disorder, for N — > 00 it 
converges to its mean valued. For the calculations below 
we have used samples of length 5 x 10 5 for relatively large 
values of W and up to 4 x 10 6 for small values of W . 



III. SCALING APPROACH FOR THE 
EIGENSTATES 

Our interest is dedicated to the structure of eigenstates 
for infinite as well as finite samples, as we tune the dis- 
order and the electric field strength. Unlike the simpler 
case of infinite samples, however, the meaning of a lo- 
calization length for finite samples is not clear. Below 
we follow the approach developed in the theory of quasi- 
1D disordered solids which is based on tha-evaluation of 
multifractal localization lengths (see, e.g.,E3). The great 
advantage of this approach is the applicability in both 
finite and infinite samples. 

One of the commonly used quantities in this ap- 
proach is the so-called entropic localization length, de- 
fined through the information entropy Hn of eigenstates, 



iV 



H N = w n ln w n ; 



1^1 



(6) 



where ip n is the n— th component of an eigenstate in a 
given basis. For eigenstates normalized as ^2 n w n = 1, 



2 



the simplest case of (p n — N^ 1 / 2 results in an entropy 
equal to the maximum value: Hn = hi (AT) . We there- 
fore define the localization length L as the number of 
basis states occupied by the eigenstate <p n ; the latter is 
equal to exp(7iAr). We note that in general, the ampli- 
tudes ip n fluctuate strongly with n and thus the coeffi- 
cient of proportionality between L and Z^ depends on 
the type of fluctuations. 

In order to study the properties of eigenstates in quasi- 
lD solids, localized on some scale in the finite basis, it 
was suggested ir£3 to normalize the localization length L 
in such a way that in the extreme case of fully extended 
states the quantity L is equal to the size of the basis N. 
In such an approach, the entropic localization length L\ 
is defined as 

Li = N exp(< H N > - H ief ) (7) 

In Eq. (0) the ensemble average < ... > is performed 
over the number of eigenstates with the same structure 
and over realizations of the disorder potential. The nor- 
malization factor 7i lc f has the meaning of an average en- 
tropy of the completely extended random eigenstates in 
the finite basis. For the quasi-ID case the distribution of 
components ip n is assumed to correspond to the Gaussian 
Orthogonal Ensemble (GOE) E3. 

Analogously, a whole set of localization lengths L q can 
be defined in the following wayE3: 



Our main question is whether relations of the type of 



L q = N ( 



where P^) is the average value of P q for the reference en- 
semble of completely extended states. One should note 
that for the particular case q = 2 the quantity Pi is 
known as the participation ratio. 

In the context of quasi- ID disordered models in the 
presence of constant electric fielcoESEJ, it was shown that 
all global properties of eigenfunctions are described by 
the following localization parameters 
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(9) 



where the superscript oo (N) denotes infinite (finite) sam- 
ples, respectively. Moreover it was found that P™' N obey 
some scaling law, i.e. they depend only on the ratio of 
the characteristic lengths of the system. In the case of 
infinite samples only two length scales, i.e. Zoo and Z e i, are 
relevant. For finite samples, however, a third length N, 
which is the actual size of the sample, comes into play 
and has to be taken into account in the scaling theory. 
According to the scaling conjecture in the modern the- 
ory of disordered solids, it was found tb|at|£0|r quasi- ID 
systems in the presence of electric fielcffij'E-il the f3° 
follow the scaling laws: 



JV 



= /3£°(Aoo); q — P q (Aoo, Ajy), 

1 ^OO \ Zoo 

where Aoo = t— ; Ajy = — • 



(10) 



Eq. (10) are also applicable for our ID AncLeESop tight- 
binding model with electric field. In RefsJliEil is was 
shown analytically that the eigenstates in ID and quasi- 
ID disordered systems without electric field, possess the 
same gross structure (envelope) on scales comparable 
with the localization length, while their statistical prop- 
erties are quite different on a much finer scale of the order 
of the lattice spacing. That is the reason why many scal- 
ing laws, which are dominated by the fluctuations of the 
envelope, hold both for ID and quasi-ID. The validity of 
such a statement is however questionable in the presence 
of electric fiels. We will show that such a similarity be- 
tween quasi-ID and ID disordered systems persist also 
in this case. 

The first nontrivial question in this context arises 
about the reference ensemble for the computation of the 
average entropy H r .f, .Indeed, in application to ID An- 
derson type modelailTEj the reference ensemble cannot 
be chosen as an ensemble of full random matrices, like 
the GOE. This point is related to the fact that in the ID 
tight-binding case fully extended states are not Gaussian 
random functions but just plane waves which arise for 
zero disorder. In the presence of electric field, the situa- 
tion is even more complicated due to strong dependence 
of the eigenstates on the electric field. However, and this 
is our expectation, in spite of the above mentioned dif- 
ference^ scaling properties of the eigenstates of the ID 
model (p]) are of the generic type discovered for quasi-ID 
systems. 

For this reason and in the spirit of Refs.ll£llI3, we define 
the normalization factors 7i re f and from the solution 
of Eq. (||) for zero disorder and electric field, i.e. e„ = 
and F = 



E k 



2V cos ■ 



kit 
N + l 



N + l 



nkn 



with k,n = 1, . . . , N. The entropy Tt lc f and the partici- 

(2) 

pation ration P^J of the above eigenstates in the large N 
limit has the same value for every eigenvalue E k , i.e. 



H Tei = ln(2AT) - 1 ; P^ = 



2N' 



(12) 



IV. SCALING PROPERTIES OF LOCALIZATION 
LENGTHS 



A. Infinite samples 



In this Section we analyze the scaling properties of 
eigenstates of infinite systems. In numerical studies the 
matrices are obviously of finite size N. However in our 



3 



analysis below we will always consider the case that 
N 3> loo, lei, an( i thus the finite (but large) size of the 
matrix becomes irrelevant. We therefore have used these 
data to investigate our scaling assumption for /3^° (see 
Eq.®). 

As was mentioned in Section II the introduction of a 
non-zero electric field F ^ 0, results in an additional 
length scale l e \. This length arises when we consider a 
cross section of the energy band locally tilted by the elec- 
tric field: — V/2 + Fn < E < V/2 + Fn for an energy level 
E. Therefore the scaling parameter Aoo = loo/l c \ enters 
the problem. Furthermore, if we consider the localization 
lengths ZA') of Eqs. (0),(|j) as the typical length, which 
contains most of an eigenvector normalization, we expect 
that 



Ac 
A. 



< 1 
> 1 



(13) 



i.e. the exponentially localized states progressively be- 
comes localized factorially as the field increases. This 
is due to the fact that, for weak electric field we have 
Aoo *C 1, and thus the dominant localization mechanism, 
i.e. the one that produces the shortest localization length 
scale, is the one related to the randomness. From now 
on we will refer to this as the " Anderson regime" . In 
the opposite limit Aoo 1, the dominant localization 
mechanism is due to the electric field. We will refer to 
this regime as the " Stark regime" . Based on the previous 
considerations we expect that the L g 's have the following 
scaling form, (see also RefJj for an equivalent reasoning 
for quasi-lD systems) 



/(Aoo) with /(Aoo) 



1 A c 



< 1 
> 1 



(14) 



where /(Aoo) is related to the scaling function (3^° as 
= Aoo/(Aoo)- 
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FIG. 2. One-parameter scaling of the (nearly) infinite sam- 
ple (N = 10 4 ) upon variation of Aoo in a range N 2> £ c i,/oo 
(W = 2.62,3.87,5, 7.35, 10, eF G [10~ 4 ,2]). A least squares fit 
according to Eq. ( p^ ) is shown as dashed line, (a) Scaling of /3J° 
(oj = 4.45 and a\ = 0.55). (b) Scaling of (ag = 4.43 and 
a\ = 0.37). 



Our aim in this paragraph is to support the above men- 
tioned scaling law based on numerical data and to extend 
our knowledge on the structure of the eigenstates in the 
intermediate regime between the two discussed limits. In 
order to study scaling properties of the localized eigen- 
states we have used the transfer matrix method for the 
calculation of Zoo as well as the direct diagonalization of 
the Hamiltonians that are associated with Eq. ([!]), for 
finite but long chains of size N = 10 4 . In all numerical 
calculations below we used V = 1. We then varied the 
disorder strength W as well as the dc field strength in a 
regime, where always Zoo, lei ^ N. One should stress here 
that both localization lengths Zoo and L q are functions of 
the energy E. For this reason, in our numerical exper- 
iments we consider ensembles of states specified by the 
values of the energy E in a small window E £ [0.95, 1.05]. 
The size of the energy window was chosen in such a way 
that the localization length Z^ is approximately constant 
inside this window (in all the cases 4^ < 0-06 ). The 

v toe 

values of pf 3 and f3%° are then obtained by averaging over 
all eigenstates which are found inside the energy window 
for a set of different realizations of disorder. As a result, 
the total number of eigenstates used for the calculation 
of (3^° were more than 1500. 

A detailed analysis of the numerical data gives evi- 
dence for a scaling behaviour of the form of Eqs. ( |l0| ) , ( |l4| ) 
with the scaling function 



/E° « a°Jl - exp(-XAoo)) 



(15) 



where the parameters a°, 
squares fit. We have found that aP q = 4.45 (4.43) and 
a q = 0.55 (0.37) for q = 1 (2). We notice here, that a 
similar scaling function was-Jjound in the framework of 
quasi- ID systems for q = lliil. Our data together with 
a fit according to Eq. ( |l5| ) are presented in Fig. 2. We 
observe a nice agreement with the scaling assumption of 
Eqs. 



B. Finite samples 



In realistic situations one always deals with finite sam- 
ples. In such cases the understanding of the statistical 
properties of conductance are of major importance. Since 
these properties are directly related to the structure of 
eigenstates, it is important to investigate the statistical 
properties of eigenstates for finite systems. This is the 
goal of the present subsection. 

..finite N and zero electric field, it was shown 
3, that the statistical properties of ID Anderson 
models are characterized by a single scaling parameter 
Aw = loo/N. Moreover, the scaling relation for the eigen- 
vectors was found to be very simple 
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c q Xn 



1 + c„A 



N 



(16) 



where the constants c q were found to be c\ ~ 2.6 and 
C2 ~ 1.5. In fact, this scaling relation is exact only for 
q — 2. For other cases of small values q, including 5=1, 
however, it is very close to the correct one (see details 
inB). 
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FIG. 3. Finite sample scaling of /3y as a function of Aoo with 
Ajv = 1. Different symbols correspond to various sample sizes JV G 
[200,1000] and disorder strengths W G [0.3,0.6]; the electric field 
was chosen appropriately (eF € [5 ■ 10 — 4 , 5 ■ 10 -1 ]). Dashed lines 
correspond to Eq. ([rij), where the values c q were taken from a least 
squares fit of Fig. 4a, b (see below). Full lines represent the scaling 
law derived in Eq. (|ig|), where the values for aP q where taken from 
Fig. 2. (a) Scaling of 0f (a° = 4.45, ci = 2.59). (b) Scaling of 
(o§ = 4.43, c 2 = 1.45). 



Once the electric field is turned on, however, a new 
length scale l c \ (with respect to the standard Anderson 
models) appears. This can be seen already from the pre- 
vious paragraph, where on the basis of numerical results 
we were able to conclude that for infinite ID Anderson 
models in the presence of an electric field the scaling 
properties of the eigenvectors are characterized by the 
parameter Aoo. Since the sample size now enters as a 
third length scale, the second scaling parameter Ajv is 
likely to show up. Thus we expect that the statistical 
properties of the eigenstates, and accordingly the 
are going to be determined by the two parameters Aoo 
and Ajv which arise due to the competition between the 
characteristic lengths l x and l c \ of the corresponding in- 
finite sample and the actual size N of the sample. In 
the rest of the Section we are going to present numeri- 
cal evidence that for finite ID Anderson models in the 
presence of an electric field, the statistical properties of 
the eigenstates are characterized by the two scaling pa- 
rameters Aoo and Ajv- To this end, we will concentrate 
on the localization measures P% (see Eq. (g)), which are 
the finite size counterparts of /3°°. 

To find the localization lengths L q for finite samples of 
size N, we have used the same approach as in the previ- 



ous subsection. The average values of L q were calculated 
by choosing only the eigenstates which had eigenvalues 
within a small energy window E G [0.95, 1.05]. Addition- 
ally, we performed an ensemble averaging over at least 
100 realizations of the disorder potential. For each N 
the total averaging thus involved several hundreds up to 
several thousands of eigenstates. In all these calculations 
the sample size was varied from N — 200 up to 1000. 



To test the scaling assumption (10) for finite systems 
we first analyze the behaviour of the localization mea- 
sures (3^2 once Ajv is fixed. This is the finite sample 
counterpart of the scaling analysis presented in the pre- 
vious subsection. In Fig. 3 we report our numerical 
results. Different symbols correspond to various sam- 
ple sizes N and disorder strengths W such that always 
Ajv = 1- The good overlap confirms the scaling depen- 
dence (3 q = (3 q (Xoo) conjectured in Eq. (|io|). 

Let us now try to gain some insight in the asymptotic 
form of the scaling law of /3^(Aoo,A/v = const). For 
Aoc *C 1, the system is in the Anderson regime, where 
[}q is given by Eq. (|l6|). The latter expression does not 
depend on Aoo , and thus we can conclude that (3 q has to 
saturate to a constant which is given by 



/3^(Aoo < l,A/v) 



1 



(17) 



where A^r = const. The dashed lines in Fig. 3 show the 
expected saturation plateau given by ( |l7j ) for Ajv = 1. 
The agreement with the numerical data is very good. 

In the opposite limit Aoo S> 1, Stark localization sets 
the dominant length scale l e \. Since we can always choose 
the strength of the electric field F such that N, Zoo 3> 1<A, 
and assuming continuity in the form of j3 q , we can ap- 
proximate the latter with the help of Eq. (|l5|). For 
Aqo ^ 1 this formula yields /3£° « oP q . Next, by changing 
variables and going from /3£° to (3 q we get 



N 



XNCLq 
Aoo 



(18) 



Displaying j3 q versus Aoo in a double logarithmic fashion, 
this yields 



ln(/3f ) » ln(X N a° q ) ~ ln(Aoo) 



(19) 



This linear behaviour ( |l9| ) is shown by solid lines in 
Fig. 3, it describes approximately the numerical data. 
Deviations are due to the fact that the approximation 
via the scaling law of Eq. ( |l5| ) actually requires not only 
l c \ -C N but also loo -C N, which is not fulfilled in our 
case. Nevertheless it gives a reasonable estimate. 

We now turn to the case, where Aoo is fixed and Ajv 
varies. Our numerical results, for Aoo = 0.01,1,20, 
corresponding to the Anderson, intermediate and Stark 
regime, respectively, are shown in Fig. 4 where now we 
refer to the new variable 



Y„ = 



1-/3, 



(20) 
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The points corresponding to the same Aoo (but different 
loo and l \) fall onto the same smooth curve with a good 
accuracy, confirming the scaling hypothesis (|io|). From 
Fig. 4 one can see that the behaviour of Y q is different 
in the two regimes Aoo 3> 1 (Aoo <C 1) where localization 
is due to the Stark (Anderson) mechanism. As a matter 
of fact, as we are increasing Aoo two asymptotic regimes 
start to build up which have the same slope and differ 
only by a constant shift. 

To understand, the behaviour of Y q (X^) as a function 
of Aoo, we first try to illuminate the limiting cases. Let 
us start with the limit Aoo <C 1. This condition, defines 
the Anderson regime, where Eq. ( |l6| ) holds and thus a 
behaviour 



hi(Y q ) = ln(Ajv) +lnc q 



(21) 



in terms of the new variable Y q is expected over the whole 
range of A at. This expectation is shown in Fig. 4 by solid 
lines. Since Aoo is small but still finite we can estimate 
the c q also from Eq. (|l5|). The limit Ajv «C Aoo < 1 
corresponds to the Anderson regime of a nearly infinite 
sample. Expanding Eq. (|l5| ) to first order yields 0? 
a g a qAoo. Substituting this expression into Eq. (| 
end up with the following term for (3 t 



N 



oN 



a q a q AN 



(22) 



Inserting (|22|) into the definition of Y q and assuming 
A at <c 1 we get 



\n(Y q ) « \n(\ N ) + ln(ay q ) 



(23) 



which implies c q 



a q a\. 



A comparison of c q — 



2.59 (1.55) and aP q a q = 2.45 (1.64) shows a very good 
agreement. Thus, we conclude that our scaling function 
( |l5| ) is consistent with Eq. ( |2l| ) as it should be in the 
above limit. 
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FIG. 4. Scaling of K = 
variation of Ajy for Aoo = 0.01,1,20 



in the finite sample upon 
Different symbols corre- 
spond to various disorder strengths W £ [0.1,6] and sample sizes 
N 6 [200,1000]. Full and dashed lines correspond to the lim- 
iting cases given by Eq. (|l(j) (with ci = 2.59, C2 = 1.45) and 
Eq. (0) (with a° = 4.45, a% = 4.35), respectively. (a,b) Scal- 
ing of Yi,l2 for Aoo = 0.01 (Anderson regime) (c,d) Scaling 
of Yi , Y-2 for Aoo = 1 (intermediate regime) (e,f) Scaling of Yi, Yi 
for Aoo = 20 (Stark regime) 

In the opposite limit of Aoo 3> 1, we have to distinguish 
between the following two cases. When Xn ^S> Aoo, the 
sample size TV sets the smallest length scale. In this limit 
the eigenstates extend over the whole sample. Then by 
continuity we expect that the behaviour of Y q (X]\[) for 
A at 3> 1 will be given by Eq. (pi"|). Our numerical data 
(see Figs. 4e,f) support this hypothesis. The second case, 
in which A at <C Aoo holds, can be viewed as the extreme 
Stark regime of an infinite sample. In that limit one ob- 
tains Eq. (|l|) again, which yields 



ln(y,)«ln(Ajv)+lnl ^- 



(24) 
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The asymptotic behaviour ( p4j ) is reported in Figs. 4e,f 
with dashed lines and agrees quite well with our numer- 
ical data. 

From the above analysis we conclude that at 



G 



Aoo =s A w (25) 

two asymptotic regimes are created due to the interplay 
between Anderson and Stark localization mechanisms. 
Although these estimations are made only on a very 
rough level, they describe our numerical findings rather 
well. 



V. SCALING OF THE DISTRIBUTION OF 
EIGENVECTORS 

As a complement to the above analysis, we show in this 
section, that the distributions of squared components of 
eigenvectors are parametrized in the same fashion by Xn 
and Aoo ■ Again we restrict ourselves to a definite energy 
window E G [0.95,1.05], where all eigenvectors corre- 
sponding to eigenvalues in that window were computed 
for several realizations of disorder. The total number of 
eigenstates were in all cases more than 1000. 

Before examining the scaling properties of the distri- 
bution let us distinguish between the various cases that 
appear due to the competition between the three charac- 
teristic length scales N, l c \, l^ (see previous section). For 
simplicity we define these regimes only by their limiting 
cases, which read as follows: 

1. loo<l e i<N & l c \<l 00 <N (infinite sample) 

2.1 00 <N<1 C \ & N<loo<l c i (Anderson regime) 

3. JV<^ el <Z 0O & ^ cl <A<Z 00 (Stark regime) 

The first pair in this list corresponds to scaling be- 
haviour of the infinite sample, since the sample size is 
always larger than the other two competing lengths. For 
this case and q = 1, 2 we have shown already in Sec- 
tion IV. A that (3^° follows a single parameter scaling with 
respect to Aoo ■ We will show here that also the whole dis- 
tribution of eigenvector components is parametrized ac- 
cording to the same scaling parameter. The first question 
which arises for the infinite sample, is the proper nor- 
malization of the squared entries of the eigenstates w n = 
|y>„| 2 . A normalization with respect to the number of 
sites N does not seem appropriate since we are inter- 
ested in the limit N — > oo. Since the w n have to scale 
with some length, however, following our previous strat- 
egy (see Eq. (^)), we introduce the variable 

r = ln(iu n Z e i) (26) 

and investigate the distribution p(r). For our calcula- 
tions we consider matrices of size N = 10 4 , while we 
lehloa *C N and varied Aoo- For each Aoo we consid- 
ered two different disorder strengths (l^ w 6, 10) and 
adjusted the dc field strength appropriately. The results 
for Aoo = 10, 1,0.1 are shown in Fig. 5a-c in a semilog- 
arithmic plot. The assumed scaling of p(r) with Aoo is 
clearly visible. 




-15 
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FIG. 5. Scaling of the entire distribution of squared eigenvec- 
tor components p(r) with Aoo in the case of nearly infinite sam- 
ples (N = 10 4 ). Two different pairs of Zoo, Z c i <C N are de- 
noted by full lines (W = 3.49, eF = 1.666,0.166,0.016) and sym- 
bols (W = 2.62, eF = 1,0.1,0.01), while keeping A^ fixed, (a) 
Aoo = 10 (b) Aoo = 1 (c) Aoo = 0.1 

The second and third pair of conditions always involve 
the sample size N . Therefore scaling according to Aoo 
and Ajv has to be taken into account. For these cases 
renormalization with respect to the sample size is mean- 
ingful; hence we define the rescaled squared entries of 
eigenstates as 

r = \n{w n N). (27) 

Before turning to the analysis of our numerical data, let 
us first qualitatively analyze the form of p(r) . In the limit 
N ^ loo, lei the system does not feel any localization due 
to disorder or electric field. Therefore all eigenstates are 
given approximately by Eq. (|Tl|). The distribution p{r) 
is then given by 

P(r)= * ■ (28) 

ny/e r {2 — e r ) 

Plotting hx{p(r)) versus r for r <gc yields a curve with 
slope 1/2, as can be verified by expanding Eq. (jg^). 
Moreover Eq. ( p8f ) shows a sharp peak around r = 0. 

In the case where Zoo -C N, Z e i the disorder sets the 
relevant length scale and the system resembles an infi- 
nite Anderson model with exponentially localized eigen- 
states w n = exp(— |n — no|/Zoo)-rJfor small r this partic- 
ular form of eigenstates leads totil 

p{r) = loo/N. (29) 

In a semilogarithmic representation this results in a 
nearly horizontal curve of height ln(?oo/-/V) for r <C 0, 
which drops rapidly for some r > 0. 
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FIG. 6. Scaling of the entire distribution of squared eigenvector 
components p(r) in the Stark regime (Aoo = 20) in the case of fi- 
nite samples (JV £ [230, 1200], W £ [0.01, 2], eF £ [10' 3 ,1]). Full 
lines and symbols denote different sample sizes (and thus different 
strengths of disorder), while Ajv is kept fixed to a chosen value, 
(a) X N = 100 (the dashed line has slope 1/2) (b) \ N = 2 (c) 
X N = 0.05 

For the second pair of conditions, i.e. the Anderson 
regime, the scaling-properties of the distribution were al- 
ready analyzed incil. A good agreement with the limiting 
Eqs. (H),© was found. 

The new and more interesting case is the pair of con- 
ditions with label 3, where the competition between N 
and l c \ is dominant. In this case A^ 3> 1 and thus the 
localization mechanism is due to the Stark effect. The re- 
sulting distribution for some representative values of Aoo 
and Ajv are shown in Fig. 6. One can clearly see that dis- 
tributions corresponding to different sample sizes N and 
disorder strengths W, but having the same A/v and Aoo 
coincide. Moreover, as we move from higher to smaller 
values of A/v the shape of the distributions changes drasti- 
cally. In the case A at ^> 1 (corresponding to N <C ^ e i) the 
eigenstates can be considered as extended with respect 
to the sample size, and thus we obtain again Eq. (|2g| ) 
(see Fig. 6a) . The peak of the distribution broadens and 
moves to the right upon an increase of Ajv as can be seen 
from Fig. 6b. At the same time, for negative values of r 
the distribution possesses long tails. This becomes more 
and more apparent as we are moving towards the Stark 
regime (Fig. 6c). In the strong field limit the eigenstates 
are essentially localized at one site i.e. w n ~ S nm . In 
this case, one has p (r ~ ln(N)) ~ 1/N, while long tails 
appear due to factorial localization. 

We conclude this Section by noticing that the scaling 
of the distributions of squared eigenvector components 
with Aoo and Ajv implies scaling of the localization pa- 
rameters ( |l0|) for arbitrary q. 



VI. SUMMARY 



We have studied a ID tight-binding model in the pres- 
ence of a constant electric field. For such a model we can 
distinguish between two regimes where localization is due 
to totally different mechanisms. The first regime is the 
Anderson regime, which is defined through the condition 
Aoo <C I. Here the localization due to disorder is the 
dominant mechanism that controls the statistical prop- 
erties of the eigenstates. In the opposite limit, Aoo ^> I, 
the localization is due to the presence of the electric field. 
This is the Stark regime. Our numerical study deals with 
the scaling properties of the eigenstates both for infinite 
and finite samples. This study was motivated by the 
remarkable scaling law thatj-Jias— been found for quasi- 
1D models with electric fieldO'ESllj. Our results indicate 
that in both infinite and finite samples with disorder and 
electric field the eigenstates have generic properties, re- 
gardless of the dimensionality of the system, provided 
that an appropriate renormalization (with respect to the 
corresponding "extended" states) is done. Thus we show 
that the similarity between ID and quasi- ID eigenstates 
should persist also for systems with electric field. 

We found that for infinite systems the statistical prop- 
erties of the eigenstates are characterized by the single 
parameter Aoo- This conclusion was based on an exten- 
sive numerical analysis of both the localization measures 
(^) and of the whole distribution of squared eigenvector 
components. Moreover for 2 we have found a simple 
scaling law (15) that describes quite nicely our numerical 
data. This expression can be used in order to find the 
strength of the applied dc electric field once is known 
for the field free model. 

Moreover, we have performed a finite length scaling 
analysis. Our numerical analysis showed that for finite 
systems the statistical properties of the eigenstates are 
characterized by two parameters, namely Aoo, as in the 
case of infinite systems, and A at. The latter parameter 
involves the actual size of the sample which enters in 
the scaling analysis as a third length scale. We found 
that the localization measures P^Li 2 show a totally dif- 
ferent asymptotic behaviour in A at — > 0, oo as we are 
increasing the parameter Aoo- Based on some analytical 
arguments we estimated that this occurs approximately 
at A/v — Aoo- This can be used as a criterion to identify 
which localization mechanism (Anderson or Stark local- 
ization) is responsible for the structure of eigenstates. 
It will be interesting to investigate if the same asymp- 
totic behaviour also holds for higher moments q. Finally, 
we studied the whole distribution of squared eigenvector 
components and showed that it also follows the same scal- 
ing behaviour with respect to the two scaling parameters 
Aoo and A at. 

The main result of our work is the fact that scaling 
properties of eigenstates of infinite systems are described 
by one parameter scaling Aoo, whereas for finite systems 
an additional parameter Ajv is also needed. In partic- 
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ular, both localization lengths, the entropy localization 
length as well as the one defined via the inverse participa- 
tion ratio, follow the universal scaling law of Eq. ( |i~0| ) af- 
ter appropriate normalization. This is in contrast to the 
standard Anderson models without electric field, where 
only one parameter (A at) is needed to describe the scaling 
properties of eigenstates. 
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